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Abstract 

The review is devoted to the integrable properties of the Generalized Kontsevich Model 
which is supposed to be an universal matrix model to describe the conformal field theories 
with c < 1 . The careful analysis of the model with arbitrary polynomial potential of order 
p + 1 is presented. In the case of monomial potential the partition function is proved to be 
a r -function of the p -reduced Kadomtsev-Petviashvili hierarchy satisfying L__ p Virasoro 
constraint. It is shown that the deformations of the "monomial" phase to "polynomial" one 
have the natural interpretation in context of so-called equivalent hierarchies. The dynamical 
transition between equivalent integrable systems is exactly along the flows of the disper- 
sionless Kadomtsev- Petviashvili hierarchy; the coefficients of the potential are shown to be 
directly related with the flat (quasiclassical) times arising m N = 2 Landau-Ginzburg topo- 
logical model. It is proved that the partition function of a generic Generalized Kontsevich 
Model can be presented as a product of "quasiclassical" factor and non-deformed partition 
function which depends only on the sum of transformed integrable flows and flat times. The 
Virasoro constraint for solution with an arbitrary potential is shown to be a standard L_ p - 
constraint of the (equivalent) p -reduced hierarchy with the times additively corrected by 
the flat coordinates. The rich structure of the model requires the implications almost all as- 
pects of the classical integrability. Therefore, the essential details of the fermionic approach 
to Kadomtsev-Petviashvili hierarchy as well as the notions of the equivalent integrable sys- 
tems and their quasiclassical analogues are collected together in parallel with step-by-step 
investigation of the suggested universal matrix model. 
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1 Introduction 



During the last years, the matrix models play an important role in the theory of 2 -dimensional 
gravity, topological models and statistical physics (see [|lj and references therein). This paper is 
devoted to the study of particular 1-matrix model in an external matrix field which is supposed to 
be "the universal" one . The structure of the model is essentially defined by the matrix integral 
of the typical form 

Z v n [M] ~ J d Xe- TrV ^ +TrXV 'W (1.1) 

where M , X are Hermitian NxN matrices and dX ~ Ilfj=i dXy . In ( |1 . 1| ) V(X) is an arbitrary 
potential (see exact formulation below). The model with V(X) = |X 3 (the Kontsevich model) 
has been derived in as a generating function of the intersection numbers on the moduli spaces, 
i.e. by purely geometrical reasons, guided by Witten's treatment of 2 -dimensional topological 
gravity . Unfortunately, the similar interpretation of more complicated model with an arbitrary 
polynomial potential is still lacking. Actually, the same model (though in somewhat implicit 
form) have appeared for the first time in [|J inspired by more "physical" arguments [||], |J. The 
advantage of the paper [3J consists of the fact that it starts from the integrable properties of 
the model from the very beginning: || gives the clear interpretation of the Kontsevich partition 
function as a concrete solution of 2 -reduced Kadomtsev-Petviashvili (KP) hierarchy, that is, the 
Korteveg-de-Vries one. This allows to generalize the original Kontsevich model immediately. In 
J7| the partition function with an arbitrary potential has been suggested as "the universal" matrix 
model under the name of Generalized Kontsevich Model (GKM) (independently, the integral 
( |1 . 1| ) with the monomial potential of finite order has been considered in papers ||- ||10|). The 



universality of GKM is based on the following facts |f7|,pT|, fL2 |: 



(i) For monomial potential it describes properly the (sophisticated) double scaling limit of any 
multimatrix model. 

(ii) The GKM partition function with a polynomial potential of order p+1 is a r -function of the 
Kadomtsev-Petviashvili hierarchy, properly reducible at the points, associated with multi- 
matrix models, to a solution of p -reduced hierarchy. Moreover, it satisfies an additional 
equation, which reduces to conventional Virasoro constraint (string equation) for multi- 
matrix models when potential is degenerates to monomial one. 

(iii) It allows the deformations of the potential associated with a given multi-matrix model to 
potentials corresponding to other models. 

(iv) GKM with arbitrary polynomial potential is directly connected with N = 2 supersymmetric 
Landau-Ginzburg theories. 



(v) The partition function (|1.1| ) with potential V(X) ~ X 2 + nlogX describes the standard 
1-matrix model before double- scaling limit. 

(vi) Adding the negative powers of X , the model gives a particular solution of Toda Lattice 
(TL)hierarchy. 

Besides, the GKM is a non-trivial (and more or less explicit) example of solution of the Kadomtsev- 
Petviashvili hierarchy corresponding to the Riemann surface of the infinite genus. As an integrable 
system with an infinite degrees of freedom, it possesses very rich structure encoding many features 
which are absent for finite-dimensional systems: the Virasoro and, more generally, W -constraints 
do not exhaust the complexity of the model. It turns out that GKM is properly designed to describe 
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the quasiclassical (dispersionless) solutions parametrized by the coefficients of the potential V(X) 
being in the same time the exact solution of the original hierarchy. This makes the study of GKM 
very promising from the point of view of the physical applications as well as in the context of 
purely mathematical aspects concerning the integrable structures. In this paper we shall deal 
with the integrable properties of ( |1 . 1[ ) only. 



To investigate the GKM model in detail, the long way to overcome is required. The point is that 
this model, being excellent example of the explicit solution of the integrable system (KP or even 
TL hierarchy), unifies many aspects of the latter. Besides well elaborated general strategy [ |13|| , 
j|nj, [TH] to describe the above hierarchies, some more subtle notions have to be implemented. 
Originally author had a temptation to dump all the details concerning the standard material 
to appendices (including the fermionic approach to r -function). After some contemplation it 
become clear that in this case the paper will contain the only Introduction as a main body with 
a lot of appendices so the structure will be the same. Therefore, by the pedagogical reasons, 
the decision arises to arrange the things as self-consistent as possible. The paper is organized as 
follows. In the first three sections we follow the approach developed in |I3|; the material here 



(except of some details) is quite standard. In Sect. 2 we give the essentials concerning the most 
important integrable system, namely, the Kadomtsev-Petviashvili hierarchy. We discuss briefly 
the pseudo-differential calculus and introduce the notions of the Baker- Akhiezer functions as well 
as the central object, the r -function, as a solution of the evolution equations. In Sect. 3 the 
fermionic approach to realization of gloo is given . It is of importance while representing the 
solutions of the KP hierarchy in the "explicit" form in terms of the fermionic correlators. In 
Sect. 4 we represent the r -function in the specific determinant form using the fermionic approach 
introduced in the previous section. Such representation is very natural from the Grassmannian 
approach to the integrable systems (TTSJ and, what is more important for the present purposes, 
this gives the simple proof of the integrability of the GKM partition function. The Generalized 
Kontsevich Model is introduced in Sect. 5. First of all, simplify the GKM partition function by 
the standard integration over the angle variables thus obtaining the integral over the eigenvalues 
Xi, x N of the matrix X . After this, we are able to write the partition function in the 
determinant form which is the starting point to investigate its integrable properties. We prove 

that, in the case of monomial potential V(X) = ^qry ■> the GKM integral is a solution of the 
p -reduced KP hierarchy. Moreover, it satisfies, in addition, the string equation. In turn, these two 
conditions fix the solution of KP hierarchy uniquely - this is exactly the GKM partition function 
with monomial potential. The case of an arbitrary polynomial potential is more complicated 
(and more rich). It requires the notion of the equivalent hierarchies |24]] which is thoroughly 
discussed in Sect. 6. We prove that the solution with polynomial potential can be generated from 
the corresponding solution of the p -reduced KP hierarchy by the action of the Virasoro group 
and is represented as another p -reduced r -function corrected by the exponential factor with 
some quadratic form. In order to investigate in detail the nature of the transformations between 
the equivalent hierarchies one more notion is required, namely, the notion of the quasiclassical 
hierarchies which is described in Sect. 7 following the approach developed in p5|-|[29||. We show 
that the quadratic form is related with the quasiclassical r -function. Moreover, we demonstrate 
that it is possible to describe the quasiclassical hierarchy directly in terms of GKM. The last section 
contains the complete description of the GKM partition function with an arbitrary polynomial 
potential. It is proved that after redefinition of times the partition function of a generic GKM 
can be presented as a product of " quasiclassical" factor and non-deformed partition function, the 
latter being the solution of the equivalent p -reduced KP hierarchy. We show how to extract 
the genuine partition function which depends only on the sum of transformed integrable flows 
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and flat (quasiclassical) times and which satisfies the standard L_ p -constraint of the (equivalent) 
p -reduced hierarchy. 



2 Kadomtsev-Petviashvili hierarchy 
2.1 KP hierarchy: Lax equations 

Let {T} = (Ti, T 2 , . . . , Tj, . . .) be the infinite set of variables. Consider the pseudo-differential 
operator (the Lax operator) 

OO Q 

L = d + J2u i+ x(T)d^; d=— (2.1) 
i=i 01 1 

where d~ l is a formal inverse to d, i.e. d~ x °<9 = dod^ 1 = 1 ; for any function /(Ti) and any 
n > 1 

d~ n o f = Y(-iy ( n + t ~ i y- od- n - 1 (2 2) 

J ' i!(n-l)! 9T| 1 ' 

Note that the Lax operator L can be written as 

L = WodoW' 1 (2.3) 

where 

OO 

W = l + Y,MT)d- 1 (2.4) 
i=i 

and the inverse W~ l can be calculated term by term using the Leibnitz rule ( |2.2| ): an easy exercise 
gives 

W~ l = 1 - Wid' 1 + (-w 2 + w 2 )d~ 2 + (-w 3 + 2u;iw 2 - wiw[ - w^d' 3 + ... (2.5) 

where ' denotes the derivative w.r.t. 7\ . Comparing ( |2.1| ) and (|2.3|) one can find the relation 
between functions {ui} and {wi} : 

m 2 = w[ 

us = —w' 2 + wiw[ (2.6) 

M 4 = — w' 3 + Wiw' 2 + w[w 2 — w\w' x — (w'x) 2 

etc. 

Let denotes the differential parts of the pseudo-differential operators L k ; for example, 

L + = d 

L\ = d 2 + 2u 2 (2.7) 
L'l = d 3 + 3u 2 d + 3(u 3 + u' 2 ) 

etc. We use also the notation L k _ to denote the purely pseudo-differential part of L k ; evidently, 
L k = Li + L k _ . 
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By definition, the dependence of functions {ui} on time variables (Ti,T 2 , . . .) is determined by 
the Lax equations 

f)T 

Wk = [L k + ,L}; k>l (2.8) 
It can be shown that this set of equations is equivalent to zero-curvature conditions 

tr-tH^ 1 (2 - 9) 

The set of equations (|2.8j ) (or, equivalently, (|2.9|) ) is called the Kadomtsev-Petviashvili (KP) 
hierarchy. Let n = 2, m = 3 in ( |2.9| ). Using the explicit expression for the differential polynomials 
L 2 + , L z + one can easily get the simplest equation of KP hierarchy - the Kadomtsev-Petviashvili 
equation: 

A( 4 ^-l2. 2 ^-^)-3^ = (2.10) 
dT^ dT 3 2 dT l dTf) &I% v ; 

2.2 Baker- Akhiezer functions 



Evolution equations for KP hierarchy ( |2.8| ) or ( |2.9| ) are the compatibility conditions of the following 
equations: 

Lty = 

(2.11) 

d T ^ = L n ,y 

J-n t 

The function \I/(T, z) which satisfies this system is called the Baker- Akhiezer function. Introduce 
the conjugation d* = —d and put 

L* = —d + (-d)- 1 o U2 + {-8)- 2 o U3 + . . . 

(2.12) 

W* = 1 + (-d)- 1 owi + (-d)- 2 ow 2 + ... 

such that L* = — (W*)' 1 odoW* . The adjoint Baker-Akhiezer function \P*(T, z) satisfies, by 
definition, the set of equations 

(2.13) 

d Tn m* = — (l™ )*\i>* 

It can be shown that solutions of the systems ( |2.12| ), ( |2 . 1 3| are represented in the form 

oo 

$>{T,z) = W(T,d)e^ z) = e^ T ' z) Y,Wi{T)z- 1 

i=o (2.14) 

**(T,z) = W*(T, d)~ l e~^ T, * ) 

where 

oo 

t{T,z) = Y,T k z k (2.15) 

k=l 

In [13|] the following fundamental theorem has been proved: 

Let ^(T, z) , ty*(T,z) be the Baker-Akhiezer functions of the KP hierarchy. There exists the 
function r(T) such that 

T (T k _ _L) 

*{T,z) = j^^ z) 

\ (2-16) 
T (T k + -^-) 

**(T,z) = =f- e ^ 

H T fc) 
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It is not hard to see that all the functions Uj(T) i > 2 can be represented in terms of r. For 
example, 

u 2 = logr 

us = 2 + ^ -9t 3 ) log r (2.17) 

«4 = \{^ - 3^9 T2 + 2^) logr - logr) 2 

etc. Substitution of the first relation to ( |2.17| ) gives the representation of the KP equation in the 
bilinear form 

1 /<9 4 r d 2 r „d 2 T\ 1 dr /<9 3 r dr 



-T 



1 2 ' V dT? 4 dT 1 ! <9T 3 + 3 <9T| ' 3 dT x \ dT? dT, 
1 / d 2 r dr \ / 9 2 r <9r 



(2.18) 



As it turns out, it is possible to rewrite all non-linear equations of KP hierarchy as an infinite 



set of bilinear equations for the r -function in more or less compact form using the Hirota 
symbols. 

One should note that it is possible to consider more general integrable system, namely the Toda 
lattice (TL) hierarchy |14[ which can be thought as a specific "gluing" of the two KP hierarchies. 



In this case the solutions depend on the two infinite sets of times, {T k } and {T^} parametrizing 
the KP parts as well on the discrete time n which mixes the KP evolutions. The r -function of TL 



hierarchy r n (T, T) also satisfies the infinite set of bilinear equations |14]; the simplest evolution 
is described by the famous Toda equation 

d 2 r n dr n dr n _ 

The main problem is to describe the generic solutions of these hierarchies. It will be done in 
section 4. 



2.3 Reduction 

The KP hierarchy is called a p -reduced one if for some natural p > 2 the operator L p has only 
differential part, i.e. 

(L p )_ = (2.20) 
In this case = L np for any n > 1 and from ( |2.11| ), ( |2.16| ) it follows that 

1 , 



T(Tfc kz k> _ n (2.21) 



dT np r(T fc ) 

From the last relation it is clear that on the level of r -function the condition of p -reduction reads 

dr{T) 



dT 



Const -r(T) n = 1,2, ... (2.22) 



Equivalently, the relations ( |2.22| ) can be taken themselves as a definition of p -reduced hierarchy. 
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3 Free field realization of gl^ 

3.1 Free fermions and vacuum states 

Let us consider the infinite set of the fermionic modes ipi , ip* , i G Z which satisfy the usual 
anticommutation relations 

{il> i ,il>l} = 6 ij , {Mj} = {^,^} = i,jeZ (3.1) 

Totally empty (true) vacuum | + oo) is determined by relations 

if>i | + oo) = , ieZ (3.2) 

Then the n-th "vacuum" state \n) is defined as follows: 

|n) = • • • | + oo) (3.3) 

thus satisfying the conditions (which themselves can be taken as definition of such state): 

ipl\n) = , k>n; i[) k \n) =0 k < n (3.4) 

Similarly, the left (dual) n-th vacuum (n| is defined by conditions 

(n\i(jl = , k < n ; (n\ipk = , k > n (3.5) 

One can select the particular state, for example, |0) and consider the normal ordering of the 
fermions with respect to this preferred vacuum. In this case the annihilation operators are ipi , i < 
and ip* , % > and, therefore, the normal ordering is defined as follows: 

= : i>ii>) : +9(-i - 1)% (3.6) 

3.2 Boson-fermion correspondence 

It is convenient to introduce the free fermionic fields 

iKz^ElM* , V^) = E^ (3-7) 

which, in turn, can be expressed in terms of the free bosonic field <p(z) 



ip(z) = q - ip log z + i -JT Z 



-k 



k&L ^ (3.8) 



[q,p] = i ; [J m , J n ] = m5 m+nfi 
according to well known formulae 

if,( z ) = : e lip{z) : = 



e* exp ^-V) x exp (- £ , ^ 
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ip*(z) = z : e 



-iip(z) . — 



(3.11) 



= ze~^ e~^ z exp (- £ ^±z k \ x exp f£ (3 ' 10) 

Note, that under the formal Hermitian conjugation 

(z) f = z- 1 ; (J fc )t = J_ fc 

(g)t = g ; ( p )t = p 

we have the involution 

, (f(2)) f ^W (3-12) 
It can be shown that vacua |n) are eigenf unctions of the operator p: 

p\n) = n\n) ; (n\p = n(n\ (3.13) 

and zero bosonic mode shifts the vacua, i.e. changes its charge 

" e im( >\n) = \n + m) 

meZ . (3.14) 

(n\e imq = (n-m\ 
Using the definition (|3.8| ) one can show that 

. e iQ vO) . . e iP<p( w ) ■ — [ z _ w \ a @ ; e ia<t>(z)+i/3<t>(w) . (3.15) 

and, therefore, 



ip(zU*(w) = — ^— : e <v(*)-i*>M 
z — w 

w 

= : if)(z)ij)*(w) : + 



(3.16) 



z — w 

The last expression being expanded near the point w ~ z enables to rewrite the bosonic field via 
the fermionic ones: 

id z y{z) = - : i/>(z)ij>*(z) : = Ik*''*" 1 (3.17) 
z feez 

or, equivalently, the bosonic currents can be represented as bilinear combination of the fermionic 
modes: 

Jk = '■ ^i+k ■■ , kez ( 3 18 ) 

Obviously, the normal ordering in ( [3.18 ) is essential only for J = p . Using ( ft. 18 ) it is easy to 
see that 

J k \n) = 

k > , neZ (3.19) 

(n\ J- k = 

One should mention that not only the bosonic currents can be expressed as bilinear combination of 
the free fermions. Actually, this is true for the whole family of gloo generators (sometimes called 
the W\ +00 -generators); for example, one can derive analogous boson- fermion corresponding for 
the Virasoro generators: 



+k ■ (3.20) 



The bosonization formulae are very useful tool to calculate different correlators containing the 
fermionic operators. 
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4 t -functions in free field representation 



In this section the solutions of KP (more generally, Toda) hierarchy are represented in the form of 
the fermionic correlators parametrized by the infinite set of continuous variables. The fermionic 
language is very convenient for the integrable systems since it enables to represent an arbitrary 
solution in the specific determinant form. This, in turn, allows to identify the GKM partition 
function with appropriate solution of the hierarchy. 



4.1 Fermionic correlators, Wick theorem and solution of KP (TL) 
hierarchy 

Let us introduce the " Hamiltonians" 

oo oo 

H(T) = £ T k J k , H(T) = J2 T k J- k (4.1) 

k=l k=l 

where {T k } and {T k } are the infinite sets of parameters (sometimes called the sets of positive and 
negative times respectively). We define the fermionic correlators ( r -functions) with the following 
parameterization by these times 

r n (T,T\g) = (n\e H ^ge- R ^\n) ee (n\g(T,T)\n) (4.2) 

where 

9 = ■ exp{ A iA*Pj} ■ (4 3) 

with \\Aij\\ e gloo ■ In the most cases we shall write r n (T,T) for brevity. We assume that the 
(infinite) matrix ||Ay|| satisfies such a requirements that the correlator (O) is well defined. As 



an example, the matrix with almost all zero entries is suitable. The wide class of the suitable 
matrices are the Jacobian ones: Ay = for \i — j\ ^> 1 . The more general conditions can be 
found in |Tj| . One should mention that the normal ordering in (|4.3|) is taken with respect the zero 



vacuum state |0) (see (|3.6|)); it is equivalent to 

Note also that every element of the type (|4.3j) rotates the fermionic modes: 



g^g 1 = R kl ^ k ; gip*g 1 = i? ifc V* (4.4) 

with some (infinite) matrix £ GLoo . As an example, the exponentials containing the Hamil- 
tonians give the transformations 

__ __ _ __ __ _ (4.5) 

e H( - T ^(z)e- H ^ = e^ T ^H{z) ; e H ^r(z)e~ H{T) = e~^ T ^r(z) 

because of commutator relations [J k ,ip(z)] = z k tp(z) , [J k ,ip*(z)] = z~ k tp*(z) (the latter are 
simple consequence of the fermionic representation ( 3.19|) ). 



The fermionic correlators introduced above have a very specific dependence on the infinite sets 
of times {T k }, {T k } . The main statement is that the correlators (|4.2j ) solve the Toda lattice 
hierarchy; in particular, as a function of the positive times {T k } these correlators are solutions 
of the KP hierarchy: each particular solution is parametrized by the given matrix ||Ay|| . It can 
be proved in full generality using the so-called bilinear identity For the local purposes it is 



10 



enough, however, to show that the simplest equations of the mentioned hierarchies are satisfied. 
It is possible to deduce them starting directly from the fermionic correlators. Note that we shall 
deal with the only KP hierarchy in what follows. Nevertheless, as an instructive example, let us 
derive 2-dimensional Toda equation which is the first equation of the Toda hierarchy. The example 
shows the natural appearance of the determinant representations in the context of the integrable 
systems; besides, the similar technique will be used below quite extensively. 



All the correlators similar to 
is applicable; as an example 



are expressed in terms of the free fields, hence, the Wick theorem 



{n\j Jll ...j> lk g(T,T)r h ---r jk \n) 
{n\g(T } T)\n) 



det 



(n\iP ia g{T,T)r jb \n) 



(n\g(T,T)\n) 



(4.6) 



a,b=l 



This key observation gives easy way to prove that the r -function ( |4.2|) satisfies the standard Toda 
equation. Indeed, using the fermionic representation ( p . 1 8[ ) of the currents Jk together with the 
definition of the vacuum states ( |3.4j ), (|3T^ 



(4.7) 



(4.t 



) one gets 

d Ti d fJn = - {n\J x e H ^ge-^ T )j^\n) = 

= - (n|'0n-lV'n^(r,T)V'„'0n-ll n ) 

Using the Wick theorem, this expression can be written in the form 

<9 Ti <9_ T n = - — {(n|'0r l -lCfl'( T )^)l ri )( ri U^( T ^)^n'0n-ll n )- 

+(n|^-i^(T,T)V:_ 1 |n)(n|^(T,T)^|n)} . 

Recalling the definitions again one can rewrite every term in the last formula in terms of the 
r -functions and their derivatives; namely, 

and, therefore, ( |4.8| ) reduces to Toda equation 

cL cL log = 



(n\ip* n g(T,T)ilj n \n) = r n+1 , 

( n \g(T, T)Vy0n-il n ) = - d T r n 

T n+l T n-l 



(4.9) 



(4.10) 



which is equivalent to ( 2.19|) . The analogous (though more involved) calculations show that r n as 
a function of the positive times Ti,T2,T% satisfies the Kadomtsev-Petviashvili equation ( |2.18 ) for 
any fixed n . Let us stress again that the complete list of bilinear equations for the r -functions 



is represented in [13, 14 



4.2 Determinant representation of r -functions 

Here we represent an arbitrary solution of the KP hierarchy in the determinant form which is 
crucial in what follows. Let us calculate the fermionic correlator (n + N\ip(fj,N) . . .ip([ii)g\n) in 
two different ways. First of all, using the definition of the vacua and applying the Wick theorem, 
the correlator can be written in the determinant form: 

(n + N\ip(/i N )...ip(ii 1 )g\n) = (n\i/)*. . ■C+JV-iV'^iv) • • .ip(m)g\n) = 

= (n\ g \n) det <"1«H-»«W)»I"> 

np 
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On the other hand, using the boson- fermion correspondence the normal ordering ( |3.15| ), and 

the formulas (|3.13|) , ( |3.14j ), ( |3.19| ) describing the action of the different operators on the vacuum 
state (N\ , one can write 



N 

(n + N\if)(/jt N ) . . . ^Mg\n) = A( A t)(n + N\ : exp {z £ 0(^)} : g\n) 



3=1 



N oo 

A(//) [| //" (n| exp { ^ T fc J k }g\n) 



(4.12) 



i=i 



fc=i 



where in the r.h.s. the r -function appears with the specific parametrization of the positive times 



1 N 

K 3=1 



(4.13) 



The parametrization ( 4 . 1 3|) has been introduced in ||16|| . We shall call such representation of times 
the Miwa parametrization (respectively, the set {//j} is called the Miwa variables). Note that 
for iV finite only first N times 7\, . . . , T/v are functionally independent. Equivalently, only first 
N equations of the KP hierarchy have a non-trivial sense (all higher equations are functionally 
dependent on the first N ones). We shall deal with such restricted hierarchy in what follows. 
Comparing the relations (|4.11|) , ( [4.12|) one arrives to the following statement. For any finite N 
the t -functions of the KP hierarchy being written in the Miwa variables (|4.13 ) can be represented 
in the determinant form 



T n (T) = (n\g\n) 
where the canonical basis vectors 

0! can V) = * 

have the following asymptotics 

i (can) i \ i—i 



det 



A(ai) 

(wK +i _i^(/x)ff|n) 
(n\g\n) 



1, 2, 



0(1 



(4.14) 



(4.15) 



(4.16) 



Moreover, the opposite statement is true. Namely, any functions r(/ii, . . . , /ijv) of the form 



r(T) 

whose basis vectors </>i(/z) , i = 



det 4>i(fij 



n = 



1, 2, ... have the asymptotics 



<M/i) = i + o 



i N 

k 3=1 



OO 



(4.17) 



(4.18) 



solve the KP hierarchy. The set satisfying the asymptotics ( |4.18|) is naturally identified 

with the projective coordinates of a point of Grassmannian |15[]. More precisely, the vectors 
{</>i (/•*)} can be transformed to the canonical ones taking the appropriate linear combinations 
(clearly, such transformation does not change the determinant in ( |4.17| )). Then, there exists 
the element ( |4.3|) of the Grassmannian such that transformed basis vectors can be written as a 
fermionic correlators ( 4 . 1 5|) (for some fixed n) and, consequently, r(/xi, . . . , //jv) have the form 
( |4.2| ) in the Miwa parametrization (|4.13|) . To summarize, any infinite set of the vectors ( 4.18Q 
describes the particular solution of KP hierarchy via determinant form ( [4.17 ) . 
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4.3 Time derivatives 



Let us find the expression of the time derivatives dr/dT k for the r -function written in the 
determinant form ( 4.14Q . As in ( |4.13| ), we assume the finite number iV of the Miwa variables. 
Hence, only first iV times T k are functionally independent and all formulas below have a sense 
for dr/dTt, . . . , dr/dT N only. From ( gOg ) 



dT k ' A(/i) 



(n + N\i/}(jj, n ) ...i/)(fjLi)J k g\n) 



III'. 

A(/i) 



N 

{(n + N\J k tfj(fi N ) . ..if)(/ii)g\n) + J^(n + N\if>(ji N ) . . . [ip(fii), J k ] ■ ■ ■ if>(jti)g\n)} 



i=l 



(4.19) 

Since the currents J k = J2jezi J ji J j+k satisfy the commutation relations [Jk,ip(fj)] = fi k ip(fi) the 
last expression can be written in the form 



dT k ' A(/i) 



n+N-l N 

(n\ipn- ■ -C+iV-i{ Yl ^j+k}^N) ■■■^i)g\n) -r n {x)J2^i (4.20) 

j=n+N—k i=l 



where, according to the definition of the vacua ( |3.5| ), the action of J k on the state (n+ N\ reduces 
to the action of the finite number of the fermionic modes with n + N — k<j<n + N — 1. This 



fact allows to represent the expression ( 4.20 ) in the compact determinant form. Indeed, since j > 
n + N — k and k < N (i.e. j > n) it is clear that (n\ipjip* +k = and the moving of the operator 
J2]=n+N-k ^Pj^Pj+k t° the left state results to appropriate shifts of the modes tp n , • • • , VCmv-i • F° r 
example, for k = 1 one gets the only correlator (n\i/;* . . . ipn+N-2' l f J n+N9' l f J (.f J 'n) ■ ■ ■ V ; (/- t i)s , l n ) an d, 
therefore, the first term in (|4.20| ) has the determinant form similar to (|4.14j ) (with the shifted 
last row (j)^ — > (/>n+i )• It is evident that for arbitrary k < N the first term in Q4.20I) can be 
represented as the sum of the shifted determinants 



(n\g\n) " 



i (can) i \ 

i (can) / \ 
<Pm+kW 

<(can)/ \ 

m +i m) 



.(can) i \ 



,(can) 
m 



1 I \ 

(can) j \ 
m+k {UN) 

/ \ 



Acan) 
An 



(4.21) 
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Hence, one arrives to the following formula: 

,(can)/ \ 



d ( det <f)\ can \fij 



dT k \ A(^) 



A(/i) 



N 

E 

m=l 



A can) i \ 



fct(Mi)-M^^ n) ^i; 

A can) 1 \ 
0m+l (Ml) 



1 (can) i 

(Pn 



(4.22) 



Introducing the formal operator which shifts the indices of the canonical basis vectors 
one can write the final answer in more compact notations: 



d I det <fii (/J*j) 



A(^) 



1 ^ 

A 7j7j E ( Bk (Vm) - /4) det ^Vi) (4-24) 



This is the first important formula we need in what follows. As an immediate application one 
can consider the translation of the notion of p -reduced KP hierarchy to the language of the 
Grassmannian. Suppose that for some natural p > 1 the quantity /i p 0^ an ^(yu) can be expanded 
in the canonical basis vectors, i.e. for any m > 1 



Writing 



//^)(//)cSpa n {0^(/.)} 



(TV) = ^i + j2 amjfI - 

3=1 



it is easy to see that for any n > 1 the following expansion holds: 



/^0<TV) = 0K>) + E«™;<A 



Tip 



^(can) 
np—j+1 



(4.25) 



(4.26) 



(4.27) 



i=i 



Due to determinant structure in ( f4.22| ) every row containing the terms 0m+n?> — /i np 0^ an ^ gives 
non-trivial contribution — « m ,„ p _ m +i0^ an ) ; 1 < m < np (provided np < N), hence, 



dr(T) 
dT, 



np 



np 

r(T)J2a 

m=l 



m,np—m+l j 



np < N 



(4.28) 



assuming that (|4.25| ) holds. In the limit N —>■ oo this is exactly the case of p -reduced KP 
hierarchy. Hence, the conditions fl4.25p and Q4.28D are equivalent [T5| . 
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4.4 Action of the Virasoro generators 

Literally the same calculation can be performed for any W -generators. Consider, for example, 
the Virasoro generators 

then, evidently, 

L k (T)r n (T) = (n\e H ^L k (J)g\n) (4.30) 



where the fermionic Virasoro generators Lfc(J) ( |3.20|) satisfy the commutations relations 



[L k (j)M»)} = fc+1 |; + ^yV>0«) = MtMM (4-31) 

Consider the subset {L_ fc (J) , > 0} . Taking into account the identity (n + N\L_ k (J) = , k > 
one gets instead of ( |4.19| ) 

T7u~ n N 

L_ fe (T)r„(T) = - ^tV E( n + Ar iV'(^iv)---[V'W^- fc (-/)]..-V'()"i)^|n) = 

^H/fJ m =l 

Y\ a~ n N 

= ~ -XTY E A -k(»m) (n + A^OM . . . V(A*i)p|n> = (4.32) 

| <7 1 77.) N 

= nkT k T n (T) — — — A-kiPm) det (p^i^j) 

^HAv m =l 

In particular, the standard r -function of the KP hierarchy r n= o(T) = r(T) satisfies the relation 

(we shall see below that the GKM partition function corresponds exactly to the choice of -vacuum 
state). 



Similarly to (|4.25|) consider the case when for some q > 1 

A^(ii)<t>t n) ^) C Spau{^™>(/x)} (4.34) 

From ( |4.33| ) it follows that the solution of the KP hierarchy is invariant w.r.t. action of the 
corresponding Virasoro generator: 

L_,(T)r(T) = (4.35) 

In the next subsection it will shown that the GKM partition function satisfies the conditions quite 
similar to ( [Egg ) and ( jOg ). 



Relations Q4.24 ) and ( [4.33| ) are the simplest examples of W -generators acting on r -functions in 



the Miwa parametrization. Using the fermionic representation it is possible to write down the 
similar expressions for the higher generators. 
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5 Generalized Kontsevich model: Preliminary investiga- 
tion 



5.1 GKM: the definition 

Recall that the standard Hermitian one-matrix model is defined as a multiple integral over n x n 
Hermitian matrix X 

Z n [t] = J e-^^'^dX (5.1) 
where the action S(X, t) depends on infinitely many coupling constants ("the times") 



oo 



S(X,t) = Y,^X k (5.2) 
fc=i 

and the measure 

n 

dX = HdX u n2^(ReX iJ )rf(ImX ii ) (5.3) 

i=l i<j 

is chosen in such a way that the following normalization condition is fixed: 

J e-^ Trx2 dX = (2vr) n2 / 2 (5.4) 



After the integration over the angle variables [21] the partition function (|5.1|) results to n -tuple 

x n 

(27T) 



integral over the eigenvalues x±, . . . , x n of the matrix X : 

n(n — 1) 

/ v •;> . , T-r - sY >• ; i f ^ 



Z n [t] = [ -^n- [ A\x) Ue- S ^da 
llfc=i K! J i=i 



(5.5) 



where 



A(x) = U(xi - X j) ( 5 6 ) 

i>j 

is the van der Monde determinant and 

n(n — 1) 



is a volume of the group SU(n) . The partition function (|5.5f) possesses a remarkable integrability 

property: as a function of times {tk} and discrete variable n (the size of the matrix) it is a 
solution of so-called Toda chain hierarchy. In particular, the function 

r n (t) = ^- Z n [t] (5.8) 

satisfies the famous Toda equation 

<9 2 logr n r n+ ir„_ 1 

= —^r (5 - 9) 



The main object we shall discuss below is quite different one-matrix integral depending on the 
external N x N Hermitian matrix M : 



f P -S(M,Y)JY 
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where the measure is the same as in ( |5.3| ) (with n substituted by N). The explicit dependence on 
the matrix M comes from the action S(M,Y) and its quadratic part S^M, Y) ; for any Taylor 
series V(Y) we set, by definition, 

S(M, Y) = Tr [V(Y + M) - V\M)Y - V(M)\ (5.11) 

such that this action does not contain the constant and linear terms in Y . The denominator 
in ( 5.1 0|) is interpreted as a natural normalization factor and is nothing but a Gaussian integral 



determined by the quadratic part of the original action: 



S 2 (M,Y) = lim \s(M,eY) (5.12) 



It is clear that the integral ( |5.10|) depends only on the eigenvalues , . . . , fj, N of the external 
matrix M . It is more reasonable, however, to use another parametrization of the partition 
function treating it as a function of the time variables T k defined by relations 

1 N 

n-rE^ ( 5 - 13 ) 

K i=i 

- these are appropriate analogues of times entering in a definition of the standard matrix model 
( |5.1| ) []. The appearance of such variables is very natural by the reasons discussed below. 



The matrix model ( 5.1 0|) is called the Generalized Kontsevich Model. The reason for this is that 
for the special choice of potential 

V(Y) = Y 3 /3 (5.14) 



the integral ( |5.10| ) becomes the partition function of original Kontsevich model 

^ (2) JdY e -V3 Trya-TVMya 



<nW\ 



J dY e ~ TrMY2 



(5.15) 



Expression (5.15) has been derived in [^] as a representation of the generating functional of inter- 
section numbers of the stable cohomology classes on the universal moduli space, i.e . it is defined 
to be a partition function of Witten's 2c? topological gravity 0. In [Q, (see also |20|, [L9f for 
alternative derivations) it was shown that as iV — ► oo Zj^ considered as a function of time 
variables ( |5.15|) , satisfies the set of Virasoro constraints 



Ll 2 )^)=0, n>-l 



(5.16) 



(5.17) 



4 2) = \ E kT k d/dT k+2n + - J2 d 2 /dT a dT b + 

^ k Odd 4 a + 6=2n 

a, 6 odd and>0 

+ j aT a bT b + -^-5 nfi - d/dT 3+2n - 

a.b odd and~>0 

Constraints ( p,16[ ) are exactly the equations imposed on the square root of the partition 

function (|5.1|) in the double-scaling limit 

1 Nevertheless, to write the explicit expression of the partition function in times {T k } requires some additional 
job. 
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5.2 GKM in the determinant form 



After the shift of the integration variable 

X = Y + M 
the numerator in (|5.10 ) can be written in the form 

f e -S(Y,M) dY _ e Tl[V(M)-MV'(M)] pJyt(M)] 

where 



F[A] = J e- TlV{x)+TlAx dX ; A = V\M) 



(5.18) 



(5.19) 



(5.20) 



Using the integration over the angular variables of the matrix X according to [p2|,p3|, one gets 



1 r N 
F[A] = (27r) 7V ( Ar - 1 )/ 2 -^- / A^He-^+^dx, 



(5.21) 



where {Aj} and {x;} are eigenvalues of the matrices A and X respectively. Therefore, the 
function F[V'(M)} in (|5.20|) can be represented as 



F[V\M)\ 



A(V'(fi)) 



det / x^'e 



j-l V(x)+V 



'^>dx\ 



N 



(5.22) 



where A{V'(fj,)) = rii>j(^' '(/•*«) ~ Vi^j)) m accordance with the definition ( |5.6|) and the unessen- 
tial constant factor is omitted. 



Proceed now to the denominator of ( 5.1U| ) 

D V N [M] = JdY e -*WV) 



(5.23) 



Making use of SU(N) -invariance of the measure dY one can easily diagonalize M in (|5.23|) . 
Of course, this does not imply any integration over angular variables and provide no factors like 
A(Y) . Then for evaluation of (|5.23| ) it remains to use the obvious rule of Gaussian integration, 



r N 

J dY e-YZ^MYM „ Y[s-}/\M) 



(5.24) 



(a constant factor is omitted again), and substitute the explicit expression for Uij(M) . If potential 
is represented clS db formal series, 



V{Y) = E ^ 



(5.25) 



k=l 



(and thus is supposed to be analytic in Y at Y = 0) , the definition ( |5.12| ) implies that 



S 2 (M,Y) = - J>J E Tr M a YM b Y 

^ k=2 \a+b=k-2 



(5.26) 



and, consequently, 



oo , , k . . k 



n =0 ^ H'j 



k=2 a+b=k-2 
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(5.27) 



Hence, 

/ e-WVdV = ftr (^)]- 1/2 (5-28) 

and substitution of ( |5.19| ), ( |5.22|) and fl5.28| ) to ( [5.1 0| ) gives the following representation of the 
GKM partition function: 



(5.29) 



Z v [M] = A( ^y ) n{ [V"(^)]- 1/2 e v ^-^ v '^\ F[V'(M)] = 

= det $r(^-)is = i 

A(a*) 

where 



5.3 Functional relations 

The vectors fl5.30|) form a linear independent infinite set. In the generic situation the basis vectors 
determining the r -function are functionally independent since they are parametrized by arbitrary 
gloo matrix. On the contrary, in GKM case the solution is parametrized , loosely speaking, by 
the vector (the coefficients of the potential V(x) ). In this sense the solution ( |5.29| ) is degenerate; 
the degeneration results to functional relations (the constraints) on the level of the basis vectors 
which, in turn, can be considered as a definition of GKM from the Grassmannian point of view 
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Consider the model, parametrized by arbitrary polynomial potential of degree p + 1 ; p > 2 : 

p+i 

k=i K 

First of all, after multiplication of ( |5.30| ) by V'(fi) , the integration by parts gives (assuming the 
vanishing boundary conditions): 



dx 

[V"{^)] l ' 2 e v ^ v '^ f {x^V'ix) - (i - l)x l - 2 }e- y ^ +xV '^dx 



(5.32) 



i.e. 



V'0m)*Y0*) = X>^ViM - (i-l)CiW ; • = 1,2,... (5.33) 

k=l 

This relation generalizes the notion of p -reduced KP hierarchy; for the monomial potential one 
gets the condition ( |4.25j ) exactly. We shall show below (Sect. |5|) that the general constraint ( |5.33| ) 
has the natural interpretation in terms of equivalent hierarchies. 
There is another type of constraint which is a generalization of ( 14.341 ) . Indeed, 

*YM = TV)] 172 e v M-" v 'M / e - v <*WMdx es A v Q*)*lM (5.34) 
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where A v (fi) is the first-order differential operator of a special form 



i a v'ifi) "' ,oj 

+ ^ 



V"(fi)dfi ^ 2[V"(fi)] 2 
Thus, we have the functional relation 

= A v (A*)*r(A*) (5-36) 
which leads to a kind of string equation similar to (4.35). To obtain the differential (w.r.t time vari- 



ables) constraint on GKM partition function resulting from ( |5.36| ), the notion of the quasiclassical 
hierarchies is required (Sect. 0). 

5.4 GKM as a solution of KP hierarchy 

We proved that the GKM partition function ( p. 10] ) is represented in the determinant form 

Z v [M] = det ggfe (5.37) 

where the vectors $Y(fJ>) are defined by ( |5.30| ). Moreover, using the steepest descent method it 
is not hard to find the following asymptotics of the GKM basis vectors: 

$rGu) = ^(l + O^- 1 )) p-oo (5.38) 



compare with ( 4.17| ), ( |4.18| ). From above consideration it follows that, being written in Miwa 



times ( 4.13Q , the partition function ( 5.37| ) solves the KP hierarchy, i.e. 

Z[T] ~ r n (T) (5.39) 
with some (yet unknown) value of the vacuum state n (see definition (|4.2|) ). 

Before proceed further, the important remark concerning the dependence on N in the formula 
( |5.37| ) deserve mentioning. The entire set {$>Y(fi)} is certainly -independent and infinite. It is 
evident that <$>Y 's are linear independent. The r.h.s. of ( |5.37| ) naturally represents the r -function 
for an infinitely large matrix M . In order to return to the case of finite N , it is enough to require 
that all eigenvalues of M , except . . . , 0jv , tend to infinity. In this sense the partition function 
ZYj[M] is independent of N ; the entire dependence on N comes from the argument M: N is 
the quantity of finite eigenvalues of M . As a simple check of consistency, let us additionally carry 
0at to infinity in ( [5.371 ), then, according to (|5.38| ), 

det^rCu,) = GuJ^ 1 ■ det N ^Y( H ) • (1 + 0(l//x J) (5.40) 

and 

AM ~ ( / . N ) iV - 1 A JV _ 1 ( / /)(l + 0(l/ / i JV )) (5.41) 

Therefore 

Z nW\ ^ Zl_ x [M\ ■ (1 + 0(1/ia v )) (5.42) 
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This is the exact statement about the iV -dependence of the GKM partition function. In this 
sense one can claim that GKM partition function is independent on iV . Therefore, we often omit 
the subscript iV in what follows. 



As the solution of the KP hierarchy, the partition function ( p.37| ) is parametrized by the coefficients 
of the polynomial V . Since the latter depends only on the finite number of parameters, the original 
matrix integral describes very particular r -function. Therefore, the question arises whether it 
possible to write down some kind of constraints which naturally select this specific solution from 
the huge set of the typical r-functions parametrized by gl^ matrix ( |4.3|) . It turns out, 

that GKM r -function satisfies the subset of Wi +00 constraints; indeed, one can find a number 
of the differential (in KP times {T^} ) operators which annihilate the function (|5.37| ). This gives 
the invariant description of the model in the spirit of ||. The problem is to describe the action of 
these operators on the r -function which is essentially written in the Miwa variables. Of course, 
due to [Q, [T7J it is well known how to reformulate all the constraints on the level of the basis 
vectors: the complete information concerning the invariant properties of the r-functions can be 
decipher from the relations similar to (|5.36|) , (|5.33|) and vice versa; this has been demonstrated 



explicitely in sections |4.3| and [4.4| . In the case of monomial potential the invariant properties of 
the basis vectors give, indeed, the complete information (see below). It is important, however, 
that relations mentioned above are not enough to describe the non-trivial evolution of the GKM 
partition function w.r.t. deformations of the potential V (say, from the monomial to arbitrary 
polynomial of the same degree). The account of such deformations results to highly involved 
mixture of the standard KP flows and so-called quasiclassical (or dispersionless) ones. In order 
to interprets the latter evolution one needs to know the action of the operators which do not 
annihilate the r -function of GKM. The non-invariant actions can not be reformulated in terms 
of the basis vectors; the explicit formulae on the level of r-functions are required. 



5.5 GKM with monomial potential, p -reduced KP hierarchy and L_ p 
constraint 

Consider the GKM partition function in the simplest case of monomial potential V(X) = ^ : 



, , e P+ 1 dXe 1 P+ 1 J det$- (u«) 
Z (P) [M] = — r J - = ^ (5.43) 

The basis vectors 

$?V) = VS^e-^T^ 1 J x i-i e-T&^dx (5.44) 
satisfy the obvious relations 

^ $«(//) = S« Ox) - (z - (5.45) 

where 

is the Kac-Schwarz operator ||. Note that up to linear term it is proportional to the Virasoro 
operator A_ p defined in (|4.33| ). 
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We have seen already that the partition function ( |5.43| ) is a r -function of KP hierarchy. Now more 
concrete statements can be made. First of all, the GKM r -function is a solution of p -reduced 
KP hierarchy. Moreover, Z^[T] is independent of times T n 



- np ■ 



«1 =0 , „ = 1,2,... (5.48) 



dT 

u± np 



In addition, the partition function (|5.43|) satisfies the L_ p constraint: 



-^- P Z^[T] + ^S1 = o (5.49) 
p oTi 



Let us give some comments concerning the relation ( |5.48| ). Due to ( |5.38|) one should note that first 
p+1 vectors ... ,&p+i(fi) has a canonical structure (|4.16|) . Therefore, for k = p the 



formula ( 4.22[) holds if one substitutes $^ , i — 1, . . . N instead of canonical GKM vectors $- can - ) 
(see more careful discussion of this point in Sect. |5.6| ). Moreover, due to ( |5.45| ) the combination 
— n p $i does not contain the vector . Hence, from ( [4.28| ) 

dZ^P) [T] 



dT p 



(5.50) 



From general KP theory one can deduce that constraint ( |5.50| ) implies all higher relations of 
the form dr [T] = Const ■ Z^[T]. Actually, it follows from the relations (|5.45|) due to 
discussion in Sect. [D| Thus, Z^[T] is, indeed, the r -function of p-reduced KP hierarchy, i.e. 
the corresponding Lax operator satisfies the constraint 

L p = {L p ) + (5.51) 

The simple proof of more strong statement ( |5.48|) , namely, the complete independence of times 
T np , n > 1 is absent, unfortunately (see [0 and, especially, JED] for details). 



To derive the constraint (|5.49 ) one needs again the canonical structure of the GKM vectors 
$^(/i), . . . , (fi) ■ It is important that because of this fact the relations (|4.24j) (with k — 1) 

and (|03D (with k = p ) can be written in terms of {$i P ' ) } . The Kac-Schwarz operator (|5.47| ) 
coincides with the formal shift operator B(p) ( |4.23| ) due to ( |5.46| ) while A_ p (n) in ( |4.33| ) is 
represented as p(A^ p \fj,) — p) . Hence, one arrives to relations 



dz (p) 

~m = am 



i 

— J2 {A {p) (p m ) - Pm) det $! P) (^) (5-52) 
L_pZ« = = -pJ_ f] (#W-/i m ) det$S p) (^) (5.53) 

^KN m =l 

thus getting the constraint (|5.48 ). Note that the latter can be written in the form 



J2Hp-k)T k T P - k + -J2(k+p) (T k+P + -4rft,i) °L = (5.54) 
ZPkTi PkTi p+l J dT k 

To conclude, the GKM r -function with monomial potential satisfies the usual L_ p -constraint 
(the integrated version of the string equation) with the shifted times 

T k -> T k H ^—5 k)P+1 (5.55) 

p+1 y ' 
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5.6 General case: V -reduction and transformation of times 



In general situation of arbitrary polynomial of degree p one gets the following matrix model: 



Z V [T] 



Tr[V(M)-MV'(M)} 
j e -s 2 (x,M)dX 



Tv[-v(x)+xv'(M)] dx 



The partition function ( 5.56 ) can be represented in the standard determinant form: 

det$r(^) 



Z V [T] 



am 



Therefore, Z V [T] is a r -function of KP hierarchy. Its basis vectors 



1/2 e V(n)-MV'Qt) 



satisfy the relations 



p+i 



E<wM-(i-i)dM 



1,2, 



fc=l 



where is the first-order differential operator 

1 d V"'(ji) 



v n {n) d^i 2\y»{n)Y 



+ fJL 



As before, these relations impose severe restrictions on the hierarchy It can be shown 
Z V [T] satisfies the generalized Virasoro constraint 



L v Z V [T] 







where 



v 



ETr 

n>l 







V"(M)M^\dT k 2 V"{^)V"{^) 



X/ \rnt 



V"{^)-V"{ N 



+ 



d 



(5.56) 



(5.57) 



(5.58) 



(5.59) 
(5.60) 

(5.61) 
^ that 

(5.62) 
(5.63) 



For monomial potential this constraint is reduced to ( 5.54| ) while, in general, it is impossible to 
write the compact expression of ( |5.63| ) in original times (|4.13|) . Nevertheless, one can construct 
the set of new times as a linear combinations of "old" ones, , in such a way that 

the operator ( |5.63| ) can be transformed to the standard one being expressed in . The way to 
find appropriate linear combinations is as follows. From ( |5.59|) one sees that GKM basis vectors 
determine the invariant point of the Grassmannian such that 

P0*)*ifa) C Span {d^Gu)} ; V(ji) = V'(ji) (5.64) 

This condition is a natural generalization of the standard p -reduction and is called the V - 
reduction. The general ideology |15| tells us that the pseudo-differential Lax operator 

L = d + U2&- 1 + u 3 d~ 2 + ... 



(5.65) 
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corresponding to this point obeys the property 

[V(L)}_ = (5.66) 

i.e. V'(L) is a differential operator of order p. Therefore, there exists the Lax operator of KP 
hierarchy 

1 = d + u 2 d~ l + u 3 d~ 2 + ... 

(5.67) 

L^> = /i\l> 

such that 

L p = T(L) (5.68) 
and, certainly, the relation between the spectral parameters of the corresponding hierarchies is 

Ji = V 1/p (fi) (5.69) 

Now it is clear that the relevant spectral parameter is jl rather then /z . Therefore, the times 
appropriate for description of V -reduced KP hierarchy should be determined by relations 

Tk = -kft 1 = -1^-^) (5.70) 

i i 

In order to find the relation between {T k } and {T^} one introduces the notion of the residue 
operation Res . For any Laurent series F(X) = J2k Fk^ k 

Res F(\)d\ = F_! (5.71) 

It is easy to see that this operation satisfies the properties 

dF(X) „ 
Res— V^dA = 
dX 

Res Fd\G = — Res Gd\F ( 5 - 72 ) 
Res Fd\G = Res F + d\G_ + Res F_d x G + 

for any two Laurent series F(X) = F + (X) + F_(X) and G(X) = G + (X) + G_(X) where F + (F_) 
are the parts of the corresponding Laurent series containing only non-negative (negative) powers 
in A . 

Using the properties of Res one finds relations: 

~ 1 00 

T fc = -V mT m ResX m - 1 V- k/p (X)dX (5.73) 

oo 

T k = Y,T m ResX- k - 1 V m/p (X)dX (5.74) 

m=k 

Let us prove now that for arbitrary polynomial potential the GKM partition function ( |5.56 ) is 
independent of time T p : 

dZ v [T(T)} 



dT p 



(5.75) 



i.e. ^'-reduced KP hierarchy resembles the standard p -reduction while considering the evolution 
along new integrable flows T k . 
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Actually, we shall derive more general formulas for the derivatives of Z w.r.t first p times 



T p . To do so, one needs to calculate the derivatives w.r.t old times. Let us apply the 



formula (|4.22|) to the GKM partition function. Immediately the problem appears. Indeed, the 



GKM vectors ( |5.58 ) have a nice integral representation, but these are not the canonical ones 
because of asymptotics ( [5.381 ). On the other hand, the formula ( J4.22 ) is valid only for canonical 
basis vectors. The compact integral representation for ^[ can \/j) is absent for GKM. Therefore, it 
is impossible to find the matrix integral representation for the derivatives d Tk Z v which is valid 
for all times k > 1 . Fortunately, this problem disappears while considering the derivatives w.r.t. 
first p times T±, ... T p . The key point is that the first p + 1 basis vectors ^(/x), . . . , 
already have a canonical form (see ( |5.38|) ). As a corollary, one can directly use (|4 . 2 2|) with the 
simple substitution $( cara ) — > & v (fx) (i.e. without any modification) for derivatives with respect 
to these times. The derivatives w.r.t. higher times do not allow such replacements. To illustrate 
this statement, one can check the "marginal" derivative d Tp Z v which contains, for example, the 
particular term (see ( |4.22|) with k = p) 



AM 



(can Vi; 



(can) 
N+p 



(5.76) 



Obviously, the first N — 1 rows in this expression can be written in terms of the GKM vectors 
( |5.58| ). The only trouble can come from the last row. But due to asymptotics ( |5.38| ) the canonical 

(can) _ ^v +p + ( ajv+p $v_ i + i ower terms) with 

$Yr i(/ii), • • • , <x n+ ®n-i(Vn) fas well as the 



vectors can be represented in the GKM basis as & N+P 



some constant and the row with entries a N+p 



N-l Kt*l) , ■■■ > ^N+p^N-l VH'N) 

rows with lower terms) does not contribute to determinant ( p.76| ). This conclusion is true for all 
other determinants resulting to d Tp Z v . Hence, the formula ( |4.22| ) with k = p remains unchanged 

if one simply substitutes $^ instead of $. can - ) 
derivatives Z v 



The same reasoning is applied, certainly, for all 
with k < p . On the contrary, the derivative d T +i Z contains the determinant 



similar to ( |5.76|) with $ 



mation $ 



(can) 
N+p+1 



(can) / ^ 



(can) 
N+p+1 



a 



(/^at) in the last row. In this case the transfor- 

+ lower terms) results to additional term proportional 
to Z v . Evidently, the higher derivatives become more and more involved while expressing them 
through non-canonical vectors 



® N+p+1 + 



N+p+1 



Due to the reasons described above only the first p derivatives have a simple integral representa- 
tions. In this case the formula ( f4.22p gives for 1 < k < p 



dZ v [T] 

on 



Tr [V(M)-MV'(M)} 
J e S 2 (X,M)dX 



J Tr [X k - M k ] e 



Tr[-V(X)+XV'(M)] 



dX 



(5.77) 



or, in compact notations, 



_d_ 

m 



logZ y [T] = (Tr X k - TrM k ) ; 1 < k < p 



(5.78) 



where 



<n*)> 



ff(X)e TT[ - V{X)+XV ' (M)] dX 

j Ty l-V(X)+XV'(M)] dx 



(5.79) 
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Indeed, the calculation of (TrX k ) does not differ in technical details from those resulting to ( |5.29| ) 
and gives exactly the r.h.s. of ( |4.22| ) with $( can ) substituted by $ y . Using the relation ( |5.74| ) 
between old and new times it is easy to find the formulas we need: 



d 



\ogZ v [T(T)] = (Tr[V k/p {X)} + -Tr[V k/p (M)}_ 



1 < k < p 



(5.80) 



Note that in the r.h.s. of ( |5.80|) is expressed through the eigenvalues of the transformed matrix 
M , i.e. M should be substituted by the solution of equation 



V(M) = M p 



The relation (|5.75|) can be readily proved now. Indeed, 



_d_ 



logZ y [T(T)] = (Tr V'(X) — Tr V'(M)) 



(5.81) 



(5.82) 



and the r.h.s. vanishes since the expression under integral is a total derivative. We proved that, 
being written in times {TJ.} , the partition function ( 5.56p has something to do with a solution 
of p -reduced KP hierarchy. Therefore, it is natural to expect that the complicated Virasoro 
constraint flOp , (fTgg) can be simplified also being represented as a differential operator w.r.t 
{Tfc} . This expectation is almost true but slightly premature: the point is that the partition 
function Z V [T(T)] is not a r function in general. Indeed, to express the partition function (|5.56| ) 
in new times ( |5.70| ) means to substitute the spectral parameters {/Xj} entering in ( |5.57| ) by the 
(formal) solution of equation (|5.69| ). Evidently, the transformation 



\i = //(l + 0(// _1 )) 



(5.83) 



destroys the structure of the van der Monde determinant and, hence, the function Z V [M(M)\ does 
not possesses the standard form. Nevertheless, the situation can be repaired: one can extract the 
genuine r -function of the p -reduced KP hierarchy from Z V [M(M)] . To describe the procedure 
we need to elaborate the notion of the equivalent hierarchies. 



6 Equivalent hierarchies 
6.1 Definition 

Consider the spectral problem L^> = where the operator L defining the KP hierarchy has a 
standard form ( [2.1| ). For any given function / 

f(p) = £ f lf S +1 / = 1 (6.1) 



2 The fermionic approach together with the above reasoning allows to write more complicated derivatives quite 
explicitely. Without proof we represent the formula 

d 2 logZ^ / / Trffe _ TrM A / TrXm _ TrM m\ \ i<k + m<p 

OTkOT m \\ /V / / 

The generalization is evident. 
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with time independent coefficients one can construct the new L -operator 



L = f(L) (6.2) 
which has the same structure as the original one. The spectral problem now is 

lm = (6.3) 

where 

n = m (6.4) 

The new operator L ( |6.2|) determines the KP hierarchy which is called the equivalent to the orig- 
inal one [23]. Introducing the differential operators Bk = (L Y , one can construct the evolution 



equations 

§- = [Bk, L] (6.5) 

which can be considered as a definition of times {T{\ . Obviously, 

B m = B* (6.6) 

The question is what is the relation between the solutions of the equivalent hierarchies determined 
by the operators L and L . First of all, one needs to establish the explicit relationship between 
{Tj} and {Tj} . The second step is to find the r -function of the " deformed" L) -hierarchy which 
corresponds to arbitrary given function r(T) of the original L -hierarchy. This gives the precise 
mapping between the equivalent hierarchies. 

6.2 Variation of the spectral parameter 

It is evident that relation ( |6.4|) can be considered as a transformation of the original spectral 
parameter [i under the action of the Virasoro generators. Let 

oo 19 11/ 

E«*uo.) - + 5 (^) - m (6.7) 



where the differential operators Ak(jj) are determined in ( |4.31| ) . The function W(n) has the 
asymptotical behavior 



W'(fi) = ^—(l + Oifi- 1 )) ix -> oo (6.8) 



where a s is a first non-zero coefficient in the sum (|6.7]) . The exponential operator exp A(fi) can 
be disentangled as 



- {^K^+^r-p^i) (...) 

where W~ x is the function inverse to W . It is convenient to introduce the function 

/(//) = W- 1 (W^) + 1) = Jl (6.10) 

which has the Laurent expansion 

f( j 2) = l 2(l + 0(fi- 1 )) /i^OO (6.11) 
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due to (|S.8D. We describe the transformation of the spectral parameter by the formula 



e w'M a* p = w-^Wifi) + l) = f{ii) (6-12) 

We have seen that the action of the operator e A ^ is expressed in terms of the function / rather 
than W . Therefore, we shall denote below the function A(fi) entering in the definition ( |B.7| ) as 
Af(fj) keeping in mind the relation ( p.lOj ) between functions /(//) and W(p) . Note that the 
relations between the coefficients and are rather complicated. 

Introduce two sets of times 

r* = f h = -W.K h (6.13) 

i i 

where p, = /(/i) . It is easy to find the relations between these times using the residue operation: 

1 oo 

T t = -E mT m Res\ m - l r\\)d\ (6.14) 

fc m=k 

oo 

T fc = ]TT m Res A- fe -7 m (A)dA ( 6 .1 5 ) 



in complete analogy with ( |5.73|) , (|5.74|) where /(A) = V l l p {\) . Note that in the operator form 



nm) = -\nK k = n^'wi r*({A*» (6.i6) 

i i 

Since 

W"(/2)d/2 = W'{ii)dn (6.17) 
the transformation (|6.15|) can be written as 

1 d 1 



T k m) = lie w<M0PiT k m) = -TErfe)) (os) 

where / _1 is an inverse function to / []. 

6.3 r -functions of the equivalent hierarchies 

Consider the correspondence between the r-functions of the equivalent hierarchies. Let r(T) be 
a solution of the original hierarchy. Using the relations Q6.15 ) one can consider it as a function 



of times {T^} , i.e. to deal with t[T(T)] . It is reasonable to assume that the latter object has 
something to do with the equivalent hierarchy determined by the operator L fl6.2|) . Actually, 
t[T(T)] is not a solution of the equivalent hierarchy. Nevertheless, being corrected by the appro- 
priate factor, this function do determines the solution we need. Namely, one can show that the 
expression 

e' Afcm?fcfm r[T(f)] = r{f) (6-19) 



3 Note that 

r\v) - W- x [W{ft-l) 



to compare with ( 6.1C ) 
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with some definite matrix A km is a r-function of the equivalent hierarchy. The most easiest way 
to prove this statement is to consider the r -functions in the determinant form ( |4.17| ). It is clear 



that the transformation of times T — > T corresponds to the transformation of the Miwa variables 
Hi = / _1 (/Ltj) . In terms of Jli the original function r[T(T)] is not the ratio of two determinants 
(since A(//)| y-in = A(/x(/x)) is not the Van der Monde determinant in terms of {/!«}) and, 
therefore, does not correspond to any r-function. Nevertheless, the r-function of the equivalent 
hierarchy can be easily extracted. Indeed, consider the identical transformation []: 

T[nf)] " nt/w))] 172 } m (6.20) 



where r(T) as a function of times ( |6.16j ) has the determinant form Q4.17 ), i.e 



f(f) = det0 - fe) (6.21) 
A(/2) 

with the basis vectors 

Uv) = [f(nmr 1/2 HKn)) ^.22) 

By the direct calculation one can show that the prefactor in the r.h.s. of ( |6.20| ) may be represented 
in the form 

TT\t>/.. I ~.\\}l/2 _ p -\A km T k T, 



WWW 2 = e^ A ^ T - (6.23) 



where 

A km = Resf k (X)d x (f m (X)) + (6.24) 
Thus, one arrives to relation ( 6.19|) . We omit brute force derivation of the identity (|6.23|) because 



the technical details are not instructive here and represent below the "physical" proof of ( |6 . 1 9| ) on 
the level of the fermionic correlators. Such approach has two advantages: first of all, it explains 
very clearly the meaning of the identical redefinition in fl6.20p -( |B.22| ) and, besides, it describes the 
explicit transformation of the point of the Grassmannian while passing to equivalent hierarchy. 

6.4 Proof of equivalence formula 

Let us consider the identity 



T(T\g) = l[e w'QH)dfii T (r\g) (6.25) 



see ( |6.18| )). The last (rather trivial) relation can be reformulated in terms of the fermionic 



correlators as follows. Using the correspondence ( |4.12|) between the fermionic correlators and the 



By /'(/^(/i)) we mean the function calculated at the point n = f 
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r -functions one can write 



(Ole^lO) = mM---1>M9\0) = ILe ""(ft) Oft (N\^Qi N ) . . ■ ^i)ff|0) _ 

W^) . . . ^)io) Di e -^fei (iV ^ ^ i)|0> ( ^ 26) 

= (A^(^)...^(/I 1 )e L /( J >g|0) (0|e H ( T >e L H J >g|0) 
(iV|^(/I i v)...^(/I 1 )e L /( J )|0) ~ (0|e^( f )e L /( J )|0) 

where 

oo 

L /( J ) = E a ^ L -fc( J ) (6-27) 
fc=i 

similarly to ( |6.7D . Thus, we proved, that the transition to the equivalent hierarchy results to the 
following identity between r -functions of the corresponding hierarchies: 

T<f\e L f {J) a) 

r{T\g) = n ;J' L J (6.28) 

i.e. one needs to redefine the times together with the appropriate change of the point of the 
Grassmannian 

9 - e L f^g = g f (6.29) 

and perform simultaneously the renormalization of the r -function. The formula (|6.28 ) coincides 
with ( |6.20| ). Indeed, the numerator is a r -function which can be written in the determinant form 
( 6.21| ) with the basis vectors 

m = ^ . = 1,2,... (6.30) 

These vectors coincide with the previously defined ones ( |6.22|) . To show this, let us move the 
Virasoro element to the left state (0| . One can discard the adjoint action of e L/ on ip* , since 

\L-k(J)M = (^y 1 (6-31) 

and, by definition, L/ contains only the Virasoro generators with k > . Therefore, for any 
i > 1 

e L / (j)^,*_ ie -L / (J) = + lower modes ( 6 .32) 

But the negative modes annihilate the left state (0| while the positive lower modes generate the 
lower basis vectors which do not contribute to the determinant det <pi{jij) ■ Therefore, without 



5 From ( f4.31| ), (J679|) - (|6 . 1 2|) it is clear that 

e L /«^) e -L/« = (^/( M )) 1/2 <K/(m)) 
i.e. the fermions are transformed as 1/2-diffcrcntials. The inverse transformation is 

During the calculations in ( |6.26 ) we are using last of these two formulas. In fact, the only things we need is that 
tp(ji) ~ e- L ^ J )V(M) e L/(J) and (N\L f = . 
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loss of generality one may assume 



(0|p|0> V ^ (0 |p|o> 

(6.33) 

where in the last equality the original basis vectors are expressed in terms of the deformed 

spectral parameter ft via the substitution ji = f~ x {p) ■ We get the basis vectors ( |6.22|) . Note 



that the appearance of the normalization factor in the definition ( 6.22| ) is very natural due to the 



fact that the basis vectors are transformed as 1/2-differentials under the action of the Virasoro 
group (see footnote || on page |30|) . 



Thus, the only problem is to calculate the trivial r -function 

r(T|e L = (0|e H ^e L ^|0> = e L ^ • 1 (6.34) 
which corresponds to the point of Grassmannian 

go = e L ' (J) (6.35) 

It is evident that this function is the exponential of quadratic combinations of T . To find the 
explicit expression let us consider its derivative with respect to the arbitrary time : 

d^r(f\e L f) = (0|e H(f) J fe e L ' (J) |0> (6.36) 



To calculate the r.h.s. of (|6.36|) one can use the following trick. The moving of the current J : 



k 



through e"k/( J ) to the right state |0) results to arising of all lower modes Ji with % < k since 

[L_j,J fc ] = -kJ k _i (6.37) 

The positive modes annihilate the right state and do not contribute to ( |6.36| ). Then move the 
negative modes through e^ f ^ back to the left. In such commutation no positive modes arise 
(due to ( |6.37| ) again). The commutation of negative modes with e H gives the linear combination 
of times due to commutation relations [J m , J n ] = m5 n+m fi ■ At last, all negative modes annihilate 
the left state (0| and the final result is r(T|e^' / ) multiplied by the linear combination of times. 
The explicit calculation is as follows. From ( |6.37| ) 



[L /5 JGu)] = d,U—J^)) (6.38) 



and, therefore, the exponentiation gives the differential operator exp ( w\/i) ^ + ( w'fo) ) 7 ) which 
can be disentangled similarly to (|6]9|). Thus, 



e L /J(/i) e - L / = J(/(/*)) (6.39) 

The inverse transformation is 



e 



- L fJ(f,)e L f = d»r\n) J{f-\^) (6.40) 



Multiplying (|6.40|) by ji k and taking the residue one gets 

j k e L f = e L fResf k {X)J{X)d\ (6.41) 
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and the action of this operator on the right state reduces to 

J fc e L '|0> = e L ^Res/ fc (A)J ( - ) (A)rfA|0) (6.42) 
where j(~)(A) denotes the linear combination of the negative current modes in the expansion 



-i 



J(A) = Y, Jk^- 1 + E = + J { ~\\) (6.43) 

fc=0 k=— oo 

Note that J(~'(X) contains only non-negative degrees of the spectral parameter A. Recall that 
we denote (_F(A)) + the part of the Laurent series F(X) containing non-negative degrees of A. 
From QQ9D 

e L /j(-)(A) = (d x f(\)j(-\f(\))) + e L f (6.44) 

One should stress that the positive modes j( + \f(X)) do not contribute to the r.h.s. of the 
last formula since d\f(X) = 1 + 0(A _1 ) and J^ + \f(X) contain only the negative degrees of A. 
Combining (|6.42 ) and (|6.44|) one gets 



J fc e L /|0) = Resf k (\)(d x f(\)j(-\f(\)))d\e L f\0) 



oo 1 (6.45) 
= E -Res/ fe (A)rf A (r(A)) + J_ m e L /|0) 



m=i m 



After substitution of ( 6.45Q to (|6.36|) and taking into account the commutation relations 



e H ^J. m e- H ^ = J. m + mf m (6.46) 

one arrives to equation 

oo 

cL r(T|e L = r(T|e L £ A km f m (6.47) 



m=l 

where 

ck 



A km = Resf k (X)d x (f m (X)) + (6.48) 
It is easy to show that A km = A mk . From ( |6.47|) the final answer is 

<0|e^^>e L /|0> = exp{- £ A km f k f m } (6.49) 

k,m=l 



and the relation (|6.28 ) between the r -functions of the equivalent hierarchies takes the form 

r[T(T)\g] = e ~\ A ^^ r(T|e L ^) (6.50) 

6.5 Equivalent GKM 

In the GKM context the equivalent hierarchies are naturally described by the function 

li = V llv {») VQa) = V\n) (6.51) 

Applying the general formula (|6.20|) one can represent the original partition function (written in 
new times T ) as follows 



1/2 



Z V [T(f )} = ^p- I] ( '-^ ) Z v [f] = c -i z v [f] (6.52) 

A W , \1>I>, J 
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where 

Ay = ResV i/p (X)d x (V j/p (X)) + (6.53) 
and the r -function of the equivalent hierarchy is described by the matrix integral 

Z v [f] = ^Q-mtf- 1 ) 1 ' 2 e Tr[y(M) - MAr/P] / e ^l-y^^ dX (6 .54) 

The r.h.s. of ( |6.54[ ) should be expressed in terms of the matrix M . Using the relation ( |6.51| ) it is 
easy to find that 

/i = - E kt k jl k ~ p (6.55) 

P k=—oo 

p+1 

Virf-vV'in) = - E t k Ji k (6.56) 

k=—co 

where 

h ee - I ^—^ResV E ^(X)dX (6.57) 

Note that the parameters t\, ... t p+ \ are independent; they are related with the coefficients of the 
potential V and can be interpreted as the times generating some integrable evolution. Indeed, 
( |6.57| ) can be considered as a set of equations which determine the coefficients of the potential as 
the functions of these additional times. The equations ( |6.57 ) naturally arise in the dispersionless 
KP hierarchy (see below). Note also that all higher positive times are zero due to polynomiality of 
V while the negative "times" , k < 0} are complicated functions of the independent positive 
times. 

Substitution od (|6.56|) to ( |6.54| ) results to the matrix integral which depends on two sets of times: 
~ z v [fA = e Er=,*-.5i j^n^r 1 ) 1 ' 2 e~ EE "^'/e 11 ™™ 1 dx\ (6.58) 

where the coefficients of V(X) are functions of quasiclassical times ti, ... t p+ i according to ( |6 . 5 7| ) . 
From (|6.53|) it follows that Ai >np = A nPti = . Moreover, t_ p = and, due to (|5.75|) , 

& = (6.59) 
dT r 



Thus, we extract the r -function of p -reduced KP hierarchy from the general matrix integral 
( |5.56| ). The last logical step to reveal the genuine integrable object hidden in ( |5.56| ) is to consider 
the part of partition function (|6.58|) without the exponential prefactor with a linear T -dependence, 
i.e. the matrix integral 

T v[f A = ^n^-y/* e -EE f dx (660) 

This is exactly the object we need. First of all, this r -function has the standard determinant 
form 

T v [f,t) = (6.61) 
J A(/x) 
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with the basis vectors 



= ^pp-i e -^" J x l - l e- V(x)+ ^ x dx (6.62) 
satisfying to p -reduction condition 

p+i 

/^f = E ^C-i - (» - Wi-i ( 6 - 63 ) 

3=1 

as well as to Virasoro-type constraint: 

Affitf = <j>f +1 (6-64) 

where 

A{Jl) = ^—-^-^—^ + -Tkt k fi k - p (6.65) 
PH P 1 a/i 2p fx p p ^ 

The partition function ( |6.60|) possesses the remarkable properties. First of all, it is a solution of 
the p -reduced KP hierarchy, i.e. 

dr v [f,t] 



dT p 



(6.66) 



this is the corollary of ( |6.58| ). Further, the relation (|6.64[) implies that ( p.6U| ) satisfies the standard 



L_ p -constraint 



L v _ v T v [T,t) = (6.67) 



1 ~ 1 00 d 

L - P = 7T E HP - k)(T k +t k )(T p _ k +t p „ k ) + - J2(k+p)(T k+p +t p+k )-^- (6.68) 
Z P k=i P k=i ol k 

where the KP times are naturally shifted by the corresponding quasiclassical ones. We shall give 
the direct proof of this statement in Sect. pTTl 



Moreover, the form of L_ p -operator ( |6.68|) give a hint that there should exist the object depending 



only on the sum T k + t k . Indeed, this is the case. Consider the product 

Z v [f,t] = r v [f,t]r {t) (6.69) 

where T (t) is a r -functions of quasiclassical p -reduced KP hierarchy. We show in Sect |8.2j that 
Z v [T, t) depends only on the sum of KP and quasiclassical times: 

/ ° A-)z v [f,t] =0, k = l,...p (6.70) 



\dT k dt k 

Of course, Z v also satisfies the constraint ( |6.671 ) 



To prove the above statements, some essentials concerning the quasiclassical hierarchies are re- 
quired. At this stage one sees that the GKM includes almost all fundamental notions of the 
integrable theory thus boiled these ingredients together. 
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7 Quasiclassical KP hierarchy 



7.1 Basis definitions 



The general treatment of the quasiclassical limit in the theory of the integrable systems can be 
found in [p5fl-[p7l and references therein. Here we outline the description of so-called quasiclassical, 



or dispersionless KP hierarchy [28] which is appropriate limit of the standard KP hierarchy (the 



careful investigation of this limit was given in p9|). Consider the quasiclassical version of the 
L -operator | 

oo 



(7.1) 



i=l 



where the functions Ui depend on the infinite set of the time variables (£1,^2, £3, • • •) = {t} and 
the evolution along these times is determined by the Lax equations 



1,2,... 



(7.2) 



where the functions ^({t}, A) are polynomials in A , and, in complete analogy with the standard 
KP theory, are defined as a non-negative parts of the corresponding degrees of the £ -operator: 



C - c 



(7.3) 



In ( |7.2| ) the Poisson bracket {., .} is the quasiclassical analog of the commutator; for any functions 
F(h,X), G(h,X) 

\FC\- &FdG 9FdG (7 A 



It is useful to introduce the additional operator 

00 00 00 



n=l 



i=l 



(7.5) 



lei 



which satisfies the equations [] 



dM 



{£\,M} 
{C,M} 



1,2,... 



(7.6) 
(7.7) 



Originally, the differential prototype of ( |7.5| ) for KP hierarchy was introduced in [^] in order to 
describe the symmetries of the evolution equations. 

In |p5 |, [ 28 1 it was proved that there exists the function ^({t}, A) whose total derivative is given 
by 



dS = J^C\dt t + Md x C 



(7.8) 



i=l 



6 In what follows we shall use the term "operator" in order to keep the resemblance with the usual KP- 
terminology; but one should perceive, certainly, that we are dealing with the functions, not with the genuine 
operators. 

7 One should stress that "the negative times" i_ 



> are the functions of the independent set {ti , i > 0} 
which are determined by the evolution equations (7.6). They have nothing to do with the actual negative times of 
the Toda lattice hierarchy. 
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and, consequently, 



(-) 



c 



+ ' 



d x S = Md x C 



(7.9) 



The function S is a direct quasiclassical analog of the logarithm of the Baker-Akhiezer function; 
the solution to (|7.9|) can be represented in the form 



S 



CO OO 1 

n=l j=l 3 



(7.10) 



7.2 Quasiclassical r -function and ^-reduction 



The notion of the quasiclassical r -function can be introduced as follows. In |28| it was proved 
that 



dh 



i+1 



dh 



j+1 -= Res CdxC J 



i,j > 1 



(7.11) 



dtj dU 

where the residue operation is defined in ( p.71| ), (|5.72| ). Therefore, there exists some function 
whose derivatives w.r.t. t{ coincide with hi+i . By definition, the quasiclassical r -function is 
defined by relations 

d log r 



H+l 



dU 



i > 1 



(7.12) 



In |29| it was shown that r -function defined above do satisfies some dispersionless variant of the 
bilinear Hirota equations, so the definition (7.12) is reasonable. 

Let us consider p -reduced quasiclassical KP hierarchy; this means that for some natural p the 
function V = C p is a polynomial in A , i.e. 



V_ = 



One can construct the "dual" function 

Q 



p p fa. 



(7.13) 



(7.14) 



which satisfy the equation 

{V,Q} = 1 (7.15) 

as a corollary of (|7.7|). In [B^], pSI the particular case of p -reduced hierarchy has been discussed, 
namely, when the function Q({t}, A) is also polynomial in A : 



o 



(7.16) 



This constraint restricts the possible solutions of p -reduced quasiclassical KP hierarchy to very 
specific subset; as Krichever has shown |ZB| the r -function satisfies the infinite set of the quasiclas- 
sical W -constraints when ( |7. 16|) holds P|. In particular, the r -function satisfies the L_i -constraint 



p+i~ 



i=l 



i=l 



dlogr 







(7.17) 



(see the proof below) . 



8 Equations (fJ^), (fTT^ ) and (\TTEi ) are the analogues of the Douglas equations (3l[] which, in turn, are equivalent 
(H to the W -constraints in the KP hierarchy. 
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7.3 Quasiclassical times and the structure of solutions 

When constraints ( I7TT3 ), (|7.16| ) are satisfied, one can construct the solutions of the hierarchy as 
follows f26fl . Note that, evidently, 

Res C~ l d x £ = S i>0 (7.18) 

Multiplying (|7.14j ) by C p ~ l ~ l d\C and taking the residue with the help (|5.71| ) it is easy to show 
that 

U = Res C^dxQ iGZ (7.19) 

From these equations for i > one can determine (at least, in principle) the coefficients of 
Q — J2iZo an d £ as the functions of times ti, t 2 , . . . , while the same equations for i < give 
then the parametrization of the "negative times" = —\h i+ i : 



1 <91ogr 

i dt 



[7.20) 



in terms of ti,t%, .... Consider the simplest situation when Q(X) is a polynomial of the first 
order. From ( 7.14|) , ( [7.16 ) it is easily seen that such condition is equivalent to switching off all 
the times with i > p + 1 : t p+2 = t p+3 = . . . = . In this case 



Q({t},\) = ^^t p+1 A + tp 



P 



(7.21) 



and equations (|7. 19| ) are reduced to 

= (g + i)t p+ i 
i(p — i) 



ResP—(\)d\ 



i < p + 1 



(7.22) 



Equations (|7.22 ) determine the coefficients of the polynomial V as the functions of first p + 1 
times ti,t 2 , • • . ,t p+ i . 

It is easy to see that the first time t\ is contained (linearly) only in A -independent term of 
V(t, A) . Therefore, 



dV 



P + l 



P+l ! 



p dti 
The Lax equations ( |7.2| ) are reduced now to the form 

P 







l,...,p) 



(7.23) 



dV 



ox (p + i)t p+1 



(7.24) 



The rest of equations ( |7.22| ) determine the functions t_j(ti, . . . ,t p+ i) , i > 1 . It is possible to 
find the explicit time dependence straightforwardly. For example, using the equation of motion 
( |7.24j ), one gets 



-J- 



ResV j/p d x V i/p = ResV i/p d x V j/p 



(7.25) 



where one uses the properties (|5.72j ) of the residue operation. In particular, the differentiation of 
t_i leads to the simple relation (since d\V x J p = dX ) 



-ResP i/? 



l p dX = 



p—i 



(7.26) 
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After integration of these equations one arrives to relation 

<91ogr 



i(jP ~ i)Ut 



dti 2(p + l)t p+1 



p—i 



(7.27) 



which is equivalent to L_i -constraint ( [7 . 1 7| ) with t 



p+2 



0. 



Note, that without loss of generality one can choose 



P 



p + l 



(7.28) 



by the rescaling of lower times and, therefore, the main equations (7.22), ( [7.24 ) acquire the 
standard form f\ 

U = 



P 



i(jp — i) 
dti ~ ~ OX 



ResV~(X)dX; i<p+l 
i = l,...,p 



(7.29) 
(7.30) 



7.4 Comparison with GKM 

The structure of the quasiclassical hierarchy has a nice interpretation in the GKM framework. 
First of all, the prepotential V'(p) = V(p.) of GKM generates the solution of the quasiclassical 
KP hierarchy subjected the constraints 

V-^) = 0; Q(jj) ~ (jl (7.31) 

The easiest way to see this is to note that the definitions (|6.57|) and ( |7.29|) are the same. Moreover, 
all the quasiclassical ingredients are naturally reproduced. Consider the first basis vector from 
the set {4>Y{n)} defined by ( |3.62j ). Neglecting the exponential prefactor one can easily see that 
the object 

V{t,n) = y/p/if- 1 J e~ vi - x)+x » P dx (7.32) 

is a Baker- Akhiezer function of the p -reduced quasiclassical KP hierarchy (recall that the 
coefficients of V are parametrized by the quasiclassical times according to Q6.571) ). It is evident 
that ^>(t,fi) has the usual asymptotic 



9{t, ai) —+ exp ( £ t kl i k ) (l + OQT 1 ] 



(7.33) 



k=l 



(7.34) 



Using equations of motion for quasiclassical KP hierarchy 

dV = _ k/p 
dt k + 

(this is the consequence of equations (|7.30|) or, equivalently the corollary of parameterization 
( |7.29| )) one can easy to show that the Baker- Akhiezer function ( |7.32| ) satisfies the usual equations 
of the p -reduced KP hierarchy: 



p(cy + ti *(t,n) = ix p ^(t,fx) 



V k + /P (d tl )*(t,v) 



(7.35) 



9 Note al so t hat the function V does not contain the term proportional to A p due to the structure of the 
£ -operator (7.1), hence, dP/dt p = . 



10 



The use of fi instead of fi should not lead to confusion. 
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where polynomials V+ p (/j) are functions of times ti, . . . t p . Hence, the function (|7.32| ) gives 
the explicit example of exact solution. On the other hand, it is important that V+ P (n) does 
not depend on t\ for k < p and, therefore, in the corresponding equations (|7.35|) we can treat 
d/dti as a formal parameter, not an operator, i.e. it is a case of quasiclassical system. Thus, we 
see that "the quasiclassical limit" can be naturally treated in the (pseudo-differential) context of 
the standard hierarchy and quasiclassical solutions are exact solutions of the full p -reduced KP 
hierarchy restricted on the "small phase space". The exact Baker- Akhiezer function (|7.32|) gives 



the explicit solution of quasiclassical evolution equations along first p flows since the standard 
relation holds: 

*(^)=exp(Et fc/ /) V \» } (7-36) 
fc=i i~{lk) 

Evaluating the Baker- Akhiezer function ( 7.32| ) by the steepest descent method, it is possible to find 



all the derivatives of the r -function entering in the r.h.s. of ( |7.36| ). To conclude, the quasiclassical 
hierarchy is determined completely by GKM integrals. 

As a consequence of the above reasoning on can see that upper pxp diagonal minor of the matrix 
Aij ( |6.53|) (which appears for the first time in the context of the equivalent hierarchies) can be 
written as the second derivative of the quasiclassical r -function. Indeed, in the p -reduced case 
the formulas (fTTTD and (|7T^) read □ 



= ResP^Pf ; ij = 1 , . . . ,p (7.37) 

and the r.h.s. is nothing but (|6.53| ). Further, the "negative" times entering to partition function 
(|6.58| ) are also represented with the help of r due to ( |7.20|) : 



dlogTo(t) , , 

kt- k = — ; k=l,...,p (7.38) 

dt k 

Before returning to the GKM r -function we need to prove some useful statement concerning the 
homogeneity of the quasiclassical r -function. 

7.5 Homogeneity property 

Lemma 7.1 The conditions ( \7. 1 imply 



S_ = (7.39) 

Proof. Recall that d\S = M.d x C . Therefore, 

Res C\d x S = Res C\Md x C = pRes C\(^-MC 1 - p )C p - 1 d x C = Res £\Qd x T (7.40) 

Since, by definition, V _ = , one gets further 

Res C + Qd x V = Res (£\Q)_d x V = Res {C\Q_)_d x V (7.41) 
On the other hand, Res C % + d x S = Res C\d x S_ . Thus, finally, 

Res £\d x S _= Res (C\Q_)_d x V Vi = l,2,... (7.42) 

n We denote now the quasiclassical r -function as To in what follows. 
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Therefore, if Q_ = then Res £ l + (d\S)_ = for any i — 1,2, ... . The last equality is equivalent 
to 

ResX i d x S_=0 (7.43) 
and, consequently, ( |7.39[ ) holds. ■ 



Lemma 7.2 The constraint ( \7. 39j) is equivalent to homogeneity condition 

°° dt 

n=l 



dt r 



t-i 



Proof. Using the explicit representation of S one gets 



(7.44) 



Res Sd x C + = Y, tn^es C n d x C\ -J2~ h j+i Res C~ j d x C\ 

j+i 3 



n=l 



(7.45) 



But for j > Res C~ j d x O + = Res C- j d x (C - D_) = Res CrH x D = iRes D'^dxC = z<^ due 
to ( [7.181 ) and, using ( |7.11| ), we have 



Res Sd x C\ = f^t n 1+1 



n=l 



dt r 



or, equivalent ly 



Res C l + d x S_ = h i+ i — V]t, 



ra=l 



hi+i 



dhi+i 
dt n 



(7.46) 



(7.47) 



and in the case S_ = one arrives to ( |7.44j ) using the identification h i+ i = —it_i . 

Note that in terms of the r -function ( [7.1 2|) the homogeneity condition (|7.44|) has the form 



L^tn—^i — = 2 log r 



n=l 



«9t, 



(7.48) 



8 Polynomial GKM: synthesis 



8.1 L_„ -constraint 



Here we represent the proof of the Virasoro constraint ( |6.67| ), fl6.68| ). 



Let 0- can ' ) be the canonical basis vectors corresponding to GKM ones, (pf , defined by (|6.62|) . 



From the general formula ( |4.22| ) one gets the expression for the derivative w.r.t. first time 7\ 



Or 



V 



m a(/2) 



, (can) , ~ 



11) ■■■ ^N-K^n) 

rN+iifii) ■■■ <P { n+1{Iin) 



N 



.V 



- T w J2fi m 

m=l 



i.l) 



40 



Now consider the action of the Virasoro generator L_ p in accordance with ( f4.33|) . The operator 
A_p(j5) being expressed via operator A{p) ( p.65| ) 



p-i 



A~ P (n) = pA(fi) - E kt kV k p -pt P - PV 



(8.2) 



k=l 



results to relation 



-L- P (T)r x 
P 



i N p-1 N 

^HA*J fc=i fe=i m =l 



(8.3) 

where the KP times are expressed through the Miwa variables fli due to (|5.70|) . In order to 
prove the analog of constraint ( |5.54j ) one should calculate the action of the operator A(fi) on the 
canonical basis vectors starting from (|6.64j ). The vectors ( |6.62| ) are not of canonical form; let 



4>uv) = ^ + oiiji 1 - 1 + . . . 



% = 1,2, . 



(see the exact expression for a, below). Now 



lican) ry 



A{n)(f>, 



(can) i(can) 

— <Pi+l 



It is evident that 



+ (a i+ i - a,i 



4^ (mi) ... <pr>(tN) 



,(can) 



+ 



+ • • • and, therefore, 
; aii — 



U) 



1.5) 



N 



E^ m )det0f a %. 



A(/2) 



and, after substitution of (|8.6| ), ( |S.1|) to ( |8.3|) , one arrives to relation 



(can) I 



i (can) i ~ \ 



+ a N+ i t 



(8.6) 



K P P dTi 



p-i 



t = - E Hp ~ H T P-ktk + (Nt p - a N+1 )r 



v 



17) 



k=l 



The last step is to calculate the coefficients ctj in the expansion (|8.4| ). Recall, that the vectors (fiY 
( |6.62| ) are related with the original ones ( |5.58|) as follows (taking into account the relation ( |6.56| )): 



{N \ V"{fi) 1 W 



1.8) 



PM 



p-i 



see the asymptotics (|5.38|) . Using the residue technique it is easy to find that 

7 p-l 



v"(fi) ^ p+i Vi 

i.e. this term does not contribute to oti . Therefore, from (|8.8|) and ( |6.55 ) 



(8.9) 



$(p) = ji + t_! + (% - i)t p w- 2 + ... 



(8.10) 
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Hence, = (i — l)t p + t_i and, consequently, 

Nt p -a N+1 = -t-i (8.11) 

where t_i is just the corresponding derivative of the quasiclassical r -function defined by ( |7.27| ); 
using the convention ( [7.28Q it reads now 

t-i = 5- X>(p - k)t k t p _ k (8.12) 
Z P k=i 



Hence, the equation (8.6) acquires the form 



1 p_1 1 °° 8\os:t v 

k){T k +t k ){T p _ k +t p _ k ) + - Y,{k+p){T k+p +t p+k ) -| = (8.13) 

and L_ p -constraint ( |6.67| ), ( |6.68| ) is proved. 

8.2 Complete description of time dependence 

Let us bring all essential facts together. 



i) We transformed the original matrix integral ( |5.56| ) in the terms of new times ( |5.74j ) as follows: 



1/2 

Z v [T(f)\ = ^® TT (^^) Z y [f] = e -i Z y [T] (8.14) 

At /0 i V/'/';' / 

(ii) The matrix depends on the quasiclassical times {t{\ related with the coefficients of 

the polynomial V(X) = V'(A) by the formula ( |7.29[ ) and is compactly written in the form 



A^t) = ResT l/p (\)d x Vi /p (\) (8.15) 

Moreover, for 1 < i, j < p 

d 2 log r (t) 

A - {t) = dUdt 3 (8 ' 16) 

and To is a r -function of the p -reduced quasiclassical KP hierarchy. 

(iii) The "preliminary" p -reduced r -function of GKM is a solution of the equivalent hierarchy; 
it is represented as 

z v \f] = ^no^r 1 ) 172 e Tr[y(M) - w ' (M)] / e ^ dx (817) 

where pP ee V'(fi) . Note that the r.h.s. of ( |S.17| ) is expressed in terms of the matrix M . In 
order to stress such dependence we denote the corresponding terms as M|~ etc. below. 



(iv) The time derivatives of the partition function ( |8.14j ) can be also written as a matrix integrals 
-L\ogZ v [T{f)] = (Tr [P k / p (X)] + - Tr [P fc/p (M)] + > 1 < k < p (8.18) 
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Let us differentiate Z v w.r.t. quasiclassical times t k keeping V'(M) = M p fixed. It is clear that 



J-Tr^(X) = -Tr[P h/p (X)]_ 
ot k 



(8.19) 



due to quasiclassical evolution equations. To calculate the derivative of {V(M) — MV'(M)}|~ 
one should take into account that, besides of the coefficient of V , the elements of the matrix M 
are also dependent on times {t k } as a functions of M . Therefore, 



9 T 



V(M)|_ -MLM P 

M M 

dV(M) DM DM 



dM dt k 



dt k 



r{% dt k 

The second term in this expression is vanished identically, thus, 

-^-Tr [V(M) — MM P )} = - Tr [r k/p (M)\ 

and, finally, 

d 



M p 



M 



dt, 



log Z V [T] = (Tr [V k/p {X)} + - Tr [V k/p (M)} + | _) 1 < k < p 



Therefore, comparing the last relation with ( p. 18 ) it follows that 



From (|8.14 ) and ( |B.23| ) one obtains 

( d d 



\ogZ v [f,t] = X> fei T, l<k<p 



\dT k dt k/ i=1 
Introduce the r -function t v in agreement with fl6.58|) , ( |6.60|) : 



T,Zi u - iTi 



where negative times t. 



Z v [T,t] = e^ i=1 r v [T,t] 
i = 1, 2, . . . satisfy the equations 
.dt-i 



A 



ki 



k = 1, ... ,p 



in accordance with (|7.25| ). Substitution of t v to (|8.24|) results to equation 



d 8 



dt. 



,dT k dt k 



^ ) log r v [f, t] = - kt. h + £ (A kl + i U -^±)f % 



t=i 



-kt. 



d log T (t) 

dt k 



due to relations ( |8.26|) . Thus, the partition function 

Z v [f,t) = r v [f,t]r (t) 



(8.20) 



(8.21) 



(8.22) 



(8.23) 



(8.24) 



(8.25) 



(8.26) 



(8.27) 



(8.28) 
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depends only on the sum T k + t k 



dZ 



dZ v 



(8.29) 
dT k dt k 

Let us find the relation between partition functions Z v and Z . Due to ( |8.14j ), ( |8.25| ) and ( 8.28| ) 
these functions are proportional up to the exponential with 



AijTiTj + it-iTi - log T (t) 



(8.30) 



Due to homogeneity relation (|7.44j ) the expression ( |3.30|) can be written as 



hence, 



-;E4W (fi+tiHfj+tj) 



l 



Z v [T(f)} = Z v (f + t)exp{--J2A 3 (t)(T t +U)(f J +t J ) 



'j 



(8.31) 



(8.32) 



This formula gives the complete description of the polynomial GKM w.r.t "quantum" ({T k }) and 
quasiclassical ({t k }) times. 
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